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Abstract 

It is found by Cai and Su that the interaction between dark energy and cold 
dark matter is likely to change the sign during the cosmological evolution. Moti- 
vated by this, we suggest a new form of interaction between dark energy and dark 
matter, which changes from negative to positive as the expansion of our universe 



(3J[) changes from decelerated to accelerated. We find that the interacting model is con- 

sistent with the second law of thermodynamics and the observational constraints. 
Tj- ■ And, we also discuss the unified adiabatic-squared sound speed of the model. 

>; 

^ ; PACS: 95.36. +x, 98.80.-k, 95.35+d 

{Sj ■ Key words: dark energy, dark matter, interaction 

Q 1 Introduction 

o 

Increasing astronomical observations [H El E] tell us that the present universe is domi- 
nated by the so-called dark energy (DE), which accounts for ~ 70% of the critical mass 
density and has been pushing the universe into accelerated expansion jU [5] . The other 
main component in the universe is cold dark matter (CDM), which accounts for ~ 30% of 
the critical mass density and behaves as the pressureless dust. However, we have known 
little about the nature of dark energy and dark matter so far. The simplest candidate 
for dark energy is Einstein's cosmological constant, which can fit the observations well so 
far. But, the cosmological constant is plagued with well-known fine-tuning and cosmic 
coincidence problems. 
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To alleviate these problems, dynamical DE models have been considered in the lit- 
erature. The simplest one is scalar-field dark energy models. So far, a wide variety of 
scalar-field dark energy models have been proposed, such as quintessence [6], phantom 
[7], /c-essence [8], tachyon [9], quintom [10], hessence [IT], etc. Other dynamical dark 
energy models include Chaplygin gas models [12], braneworld models [13], holographic 
models [JJ], agegraphic models [T5], etc. A lot of efforts have been made to explore the 
nature of dark energy. Furthermore, since no known symmetry in nature prevents or 
suppresses a nonminimal coupling between DE and CDM, there may exist interactions 
between the two components. At the same time, from the observation side, no piece of 
evidence has been so far presented against such interactions. Indeed, possible interac- 
tions between the two dark components have been discussed intensively in recent years. 
It is found that a suitable interaction can help to alleviate the coincidence problem [16] . 
Different interacting models of dark energy have been investigated [T7| [TBI . 

In the literature, the model with interaction between DE and CDM is usually de- 
scribed by the Friedmann equation 

# 2 = y(p m + p d ), K 2 = 8nG. (1) 
and the two conservation laws 

p m + 3Hp m = Q, (2) 

p d + 3H(p d + Pd ) = -Q, (3) 

where Q denotes the phenomenological interaction term; p m and p d are the energy den- 
sities of CDM and DE respectively; p d is the pressure density of DE; H = a/a is the 
Hubble parameter; a is the scale factor in the Friedmann-Robertson- Walker (FRW) met- 
ric; a dot denotes the derivative with respect to the cosmic time t. Usually, three forms 
of Q are used 

Q 1 = 3bH Pd , (4) 
Q 2 = 3bH{p d + p m ) 1 (5) 
Q 3 = 3bH Pmi (6) 

where h is the coupling constant. Then positive b means that DE decays into CDM, 
while negative b means CDM decays into DE. In the cases of Q = Q\ and Q = Q2, 
negative b would lead p m to be negative in the far future. For negative b in the case of 
Q = Q3, no such difficulty exists. But in Ref.[19j. from the thermodynamical view, it 
is argued that the second law of thermodynamics strongly favors that DE decays into 
CDM, i.e. b is positive (see Ref.[20] for a different view). So generally b is taken to be 
positive. 

However, recently it was found that the observations may favor the decaying of CDM 
into DE [20, ED [22] • Particularly, in Ref . |24] . in a way independent of specific interacting 
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forms the authors fitted the interaction term Q with observations. They found that Q 
was likely to cross the noninteracting line (Q = 0), namely the sign of interaction Q 
changed, around z = 0.5. This raises a remarkable challenge to the interacting models, 
since the usual phenomenological forms of interaction, as shown in the last paragraph, do 
not change their signs during the cosmological evolution. As noted in [23] , more general 
forms of interaction should be considered. 

In our paper, we are interested in proposing such a new form of interaction. It is 
known that our universe changes from deceleration to acceleration around z = 0.5 [25J. 
Thus the interacting term proportional to (pd — p m ) would change its sign naturally 
around z = 0.5. So we may assume the new form of interaction to be 



where a is the coupling constant and a is a positive constant of order unity. For sim- 
plicity, we take a = 1. So the new interaction term is assumed to be 



Henceforth, we will denote the new interacting model as a. The parameter a is assumed 
to be positive, since negative a would lead to negative p m in the far future. Obviously, in 
the early stage, Q4 is negative, since p m > pd- As the expansion of our universe changes 
from decelerated to accelerated, Q4 changes from negative to positive. 

Below, we first show that the interacting model with Q = Q4 is consistent with the 
second law of thermodynamics by following the argument in [T5]. Second, we calculate 
the unified adiabatic sound speeds of the interacting model. Then we compare the model 
with observations. Finally, conclusions are given. 

2 New Interaction and The Second Law of Thermo- 
dynamics 

We recall the thermodynamical description of DE and CDM in [T9J [20j [23] . A perfect 
fluid is characterized by (n, p,p, s,u a ), where n is the particle number density, p is the 
energy density, p is the pressure density, s is the entropy per particle and u a is the 
4-velocity. In the FRW universe, we take u a = (d/dt) a and u a . a = 3H. In this paper, 
a semicolon denotes the covariant derivative compatible with the FRW metric. The 
energy-momentum tensor may be assumed to be [23J 





(7) 







(9) 



Then the conservation law of the fluid u a T' 



■b=0 gives us [20] 



p + 3H(p + p) 



3ffll, 



(10) 
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where —3HU represents the phenomenological interaction between the fluid and others. 
The balance equation for the particle number is assumed to be [20] 



n + 3Hn = nT } (11) 

where T is the rate of the change of the particle number of the fluid. The temperature 
T of the fluid is defined via the Gibbs equation 

TYfa^) +*(!), (12) 



so that the variation rate of the entropy per particle is [20] 

i = ™ _ e+l r . ( i3) 

Tn Tn y ' 

By defining the entropy flow vector as 

S a = snu a , (14) 

we have [20J 

S % = ( s-^)nT-^. (15) 

(16) 



The entropy per particle is 

= p+p _ a 

8 nT T'' 
where \x is the chemical potential. Then, we have 



S% = ~^nV - 3 ™ (17) 



T T 

In the general case, the Gibbs equation can be rewritten as 

Ln l In 

and the integrability condition 

d 2 s d 2 s 



dpdn dndp 

tells us 



dn dp dp 

Using Eqs. (fTU]l . (ITT]), (jI5|) and ([15 )1 . we have 



dT dT dv 

n^- + (p + p)^- = T^f. (19) 



f „rrd P .dT dp 

f = - 3H dp- +n % +T dp- (20) 
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Generally, we should assume that at any event in spacetime, the thermodynamics-state 
of the fluid is close to the fictitious equilibrium-state at that event [23] . This implies that 
the right-hand side of Eq. (120]) is dominated by the first term, and we have approximately 

f dp 

f - - 3H i < 21 > 

Now we apply the results above to the model with interaction between DE and CDM. 
Usually the energy-momentum tensor (EMT) of CDM is taken as 

T£ = p m u\u\ + p m (g ab + u a x u\). (22) 

Then the conservation law should be 

T^, b =W b + ulu\)], h . (23) 

Here u\ is the four-velocity of CDM. In this paper, we use the subscripts m and d to 
denote the corresponding parameters of CDM and DE, respectively. With p m = and 
choosing 

U = a(p d -p m ), (24) 

we can deduce Eq.([2]) with Q = by contracting Eq. (l23|) with u\ a . Then equivalently 
we can define the effective EMT of CDM as 

= p m u\u\ + (p m - U)(g ab + u a lU \). (25) 

Obviously, the effective EMT of CDM is conserved 

rpab ri 

me;6 u# 

Similarly, although the EMT of DE 

Tf = Pd u a 2 u b 2 +p d (g ab + u a 2 u b 2 ), 

is not conserved 

Tf. b = -[U(g ab + u^] ;b , 
we can define the effective EMT of DE as 

T£ = Pd u a 2 u h 2 + (p d + U)(g ab + u%u b ), (26) 

which is also conserved 

T£. b = o. 

Here u 2 is the four- velocity of DE. We can recover Eq.Q with Q = from u 2a T£ b . b = 
if n is chosen as given in Eq. 
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The equation of state of DE is 

Pd = wp d . (27) 

In this paper, we only consider the model with constant w. For CDM, approximately 
we have [23J 

3 

Pm = n m M + -n m T m , p m = n m T m (k B = 1), (28) 

so long as T m <C M. From Eq. (l2ip . approximately we have 

T m ex a" 2 , T d oc a~ 3w . (29) 

The results tell us that as the universe expands, the temperature of CDM, T m , decreases 
and the temperature of DE, T d , increases. Thus, one may expect that at the present and 
in the future T m < T d , while in the past T m > T d . Following Ref.[T9], we assume that 
both DE and CDM have null-chemical potentials. Thus from Eq. ffT7"|) . we have 

Sl, a + S a d , a = (-L--L)Q, (30) 

The existence of interaction means that there is a transfer of energy between DE and 
CDM. It is a natural conclusion that nowadays the energy is transferred from DE to 
CDM, i.e. Q > 0, since currently T d > T m . In addition, from Eq.(13"0"]). the second law of 
thermodynamics + S d . a > and because T m < T d indicate that currently Q > 0. 

The interaction term Q = Q\ is used in the analysis in |19] : Q\ cannot change its 
sign, and is still positive even when T m > T d . It seems that at earlier, the second law 
of thermodynamics was violated, since T m > T d , Qi > and Eq.f l30p indicate that 
+ 5*^^ < and the energy is being transferred from the lower temperature DE to 
the higher temperature CDM. To overcome the difficulty, the author in [19J argued that 
the thermodynamical description breaks-down at some point, both when a < 1 and 
when a ^> 1. 

Now, we apply the analysis in [19] to the case of Q = Q4 in Eq.([8]). Obviously, 
nowadays the second law of thermodynamics is satisfied since Q4 is positive. Earlier, 
since p m > p d , Q 4 was negative, which indicates the energy is transferred from CDM to 
DE. In fact, this is just what is expected from the second law of thermodynamics when 
T m > T d . So, the interacting model a with Q = Q± is always consistent with the second 
law of thermodynamics even at early times. 

3 Adiabatic Sound Speed 

The squared sound speed cf, defined as 

s 5p' 



6 



is an important quantity for the cosmological evolution, which determines the stability 
of the cosmological evolution [26] . The adiabatic-squared sound speed c 2 a is defined as 

cl = 1 (31) 

In the interacting model a, we rewrite Eqs.([2]) and ([3]), respectively as 

p m + 3H( Pm +p^) = 0, (32) 
p d + 3H(p d + pf) = 0, (33) 

where 



p%=-<r(j> d -p m ), (34) 
Pf = Pd + cr(Pd- Pm)- (35) 



Naively, we may still define the squared sound speed of DE c 2 d as 



"sd r 

dp. 



Spd 
d 



(36) 



Actually, this is not the physical sound speed of DE. (The two physical sound speeds 
A± are shown in the Appendix.) However, we find that in the model a the stabilities of 
CDM and DE under perturbations are still determined by the squared sound speed of 
DE c 2 sd (see the Appendix for details). The corresponding adiabatic sound speed is 

eld = P ~- (37) 
Pd 

Then, using Eqs. (!27|) . we have 

eld = w. (38) 

Negative w indicates that adiabatic instabilities exist. Yet, this is not very astonishing 
since it is well-known that even in the noninteracting model with constant w, the adia- 
batic squared sound speed of DE c 2 ad is negative (c 2 ad = w < 0), leading to the adiabatic 
instabilities. 



4 Comparison with Observational Data 

In this section, first we will explore whether the interacting model a is consistent with 
the results in Ref. |27j. Second, we will explore whether the model a is consistent with 
the observational constraint on the position of the first peak of the cosmic microwave 
background power spectrum. 
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In order to explore whether the model a is consistent with the results in [27J, we 
should calculate the dimensionless coordinate distance, y(z) = H Q a r, and the two first- 
derivatives with respect to redshift, and then compare the results with the observational 
data of supernovae-type la (SN la) and radio-galaxies between the redshift z = and 
z = 1.8. Here z = — — 1 is the cosmo logical redshift and f is the radial coordinate in 
the FRW metric. 

Since dt = —a(t)dr for photons flying from their sources to observer in the flat-FRW 
universe, we have 

/ s H(z) 1 , , 

where y' = dy/dz. We define 

Pm = Anl)3i(z), Pd = PdO%{z)- (40) 

Then using Eqs.flU), (]2l), (13]), (155]) and (110]) . we have 

y ' {z) = in m0 z 1 + n d0 %Y/i> (41) 

"( \ - 3 V' T ^O^ 1 + ( 1 + U))^dO% ] ( A r )\ 

Z 1 -+- Z L ii m0 ^i -t- Sid0^2 J 

where j/" = d 2 y/dz 2 . 

Now let us calculate Si and 3 2 . By using Eqs.flS]), (]4"0j) and f[2"Tj) . the conservation 
laws 02]) and ([3]) read respectively 

rfSi 3(l + ff) 3a 

-3 1 ^1 = ^2, (43) 

da a aro 

d% 3(l + w + a) 3r (r 

— 1 is 2 = bi, (44) 

da a a 

where r = By solving the two equations, we have 



= c(l + ^ + + (1 - c)(l + z) s -, (45) 
3 2 (*) = -^-[cs+(l + z) s+ + (1 - c)S_(l + z) s -], (46) 

3(7 



where 



s± = -[2(1 + a) + w ± v^a 2 + w 2 ], (47) 

s± = ^(w± V4cr 2 + w 2 ), (48) 

1 w + ^ 

c = i g , (49) 

2 2Vw 2 + 4a 2 
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Fig. 1: y(z) versus z in the interacting model a for fixed f2 m0 = 0.272, a = 10 3 and 
different w. For comparison, the prediction of the ACDM model is also shown. 
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Fig. 2: y'(z) versus z in the interacting model a for fixed fl m0 = 0.272, a = 1CT 3 and 
different w. For comparison, the prediction of the ACDM model is also shown. 
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Fig. 3: y"(z) versus z in the interacting model a for fixed Q m o = 0.272, a = 10 3 and 
different w. For comparison, the prediction of the ACDM Model is shown, too. 

In the absence of interaction, the Eqs. (j4"Tj) and fj4"2]) reduce to Eqs.(l) and (2) of Ref. [27J. 

The evolutions of y, y' and y" with respect to z are depicted in FigJH Figj2]and FigJSJ 
We have used the different values of the equation of state parameter w, and fixed the 
values of fi m0 = 0.272 [22] and a = 10" 3 . Comparisons of Figsffl Figfjand Fig|3]with 
the corresponding figures in Ref.[2"7j reveal that the interacting model a is consistent 
with the analysis of Ref. [27] . 

Now, let us calculate the shift parameter R, which characterizes the position of the 
first peak of the cosmic microwave background spectrum and is defined as [28] 




Here z* is the redshift of decoupling. The 7-year Wilkinson Microwave Anisotropy 
Probe(WMAP) observations tell us that z* = 1091.3 ± 0.91 at la confidence level |22]. 
In this paper, we fix = 1091. The values of the shift parameter R for different sets 
of parameters in the model a and the ACDM model are displayed in Table [H We have 
fixed f2 m o = 0.272 and used the flat FRW metric. We first consider that the universe 
is filled with DE (Qdo — 0.728) and CDM (f2 m o = 0.272). Second, we consider the case 
that the universe is filled with DE (fi d0 = 0.728), CDM (fi dm0 = 0.2264), and baryon 
matter (Qbo = 0.0456) [29j. The total fractional-energy density of CDM and baryon 
matter is still fixed to be Q m0 = 0.272. The component of baryon matter is assumed to 
be evolving separately. In the ACDM model, the shift parameter R is determined by the 
total fractional energy density of CDM and baryons. The 7-year WMAP observations 
tell us R = 1.725 ± 0.018 at la confidence level. Then from the results displayed in 
Table [H we know that the interacting model a is consistent with the 7-year WMAP 
observations. 





(50) 
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Model 




R 


a 


0.272 -0.9 10~ 3 


1.712 


a 


0.272 -1.05 10" 3 


1.738 


a 


0.272 0.0456 -0.9 10~ 3 


1.713 


a 


0.272 0.0456 -1.05 10~ 3 


1.739 


ACDM 


0.272 -1 lO" 3 


1.733 



TABLE 1: The values of the shift parameter for different sets of parameters in the model 
a and the ACDM model. 

5 Conclusions 

The dark energy models with interaction between DE and CDM have been investigated 
intensively. It is argued in [SU] that the interacting models of dark energy may be key 
to solving the cosmic coincidence problem. Recently, it was found that the interaction is 
likely to cross the noninteracting line [23] . However, the usual forms of interaction used 
in the literature can not change their sign. The result in [23] raises a challenge to the 
interacting models of dark energy; more general interaction is needed to be considered. In 
[3T] , the interacting models with interaction proportional to the deceleration parameter 
are suggested. 

In this paper, we suggest the new interacting model a with interaction Q = Q4. 
Obviously the sign of Q4 changes from negative to positive as the expansion of our 
universe changes from decelerated to accelerated. We found that the interacting model 
is consistent with the second law of thermodynamics. Then we also found the squared 
sound speed of DE, c^ d , to be crucial in determining the stabilities of DE and CDM in 
the model. As in the noninteracting model with constant w, there also exist adiabatic 
instabilities in the model a with constant w due to c ad = w < 0. Furthermore, we 
compared the interacting model with observational data. And we found the interaction 
model to be consistent with the result in [27J and the 7-year WMAP observations [29J. 
Thus, we believe that the interacting model with Q = Q± is consistent with observational 
constraints. In the future, we plan to explore how to obtain the phenomenological 
interaction term Q4 from an action principle. 
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Appendix 



Here, we try to show that the squared sound speed of DE determines the stabilities of 
the cosmological evolution in the interacting model. We adopt the conformal Newtonian 
gauge. Then the perturbed metric about a sptially flat, homogeneous and isotropic FRW 
background is given to be [32J 



ds 2 = a 2 [-(l + 2<p)dT 2 + (1 - 2^)dx 2 ]., (51) 

where a is the scale factor and the perturbations of the metric are characterized by two 
potentials, <p and ip. The total EMT is 

T t ab = T m b e + T£. (52) 

The T m b e and T d a e 6 are the effective EMT of CDM and DE, respectively, 

t = p m «X+^(«X+/), (53) 
T& = PdU a A + Pf^>2 + 9 ab )- (54) 

pab _ rpab 
me;fe me;6 



Here and p'f are defined in Eqs. (l3~21) and (135]) respectively, and T^ e . b = T^ e . b = 0. 



To the first-order of perturbation, we have 



p m (r,x) = p m6 (r)[l + 5i(r,x)], (55) 
p e m ff (r, x) = pS(r) + tfpff(r, x), (56) 
<(r,x) = ^[(l - (j))(d/dT) a + d i v 1 (T^){d/dx i ) a l (57) 



and 



p d (r,x)=p ca (r)[l + 52(r,x)], (58) 
pf(r,x)=^(r) + ^ ff (r,x), (59) 
<(t,x) =a- 1 [(l-0)(9/9r) a + ^ 2 (r,x)(9/9x i ) a ]. (60) 

Here we use the subscript 6 to denote the spatially homogeneous-background value of 
the corresponding quantity, and 8\ = and 52 = ^ are the fractional perturbations 
in the energy densities of CDM and DE, respectively; v\ and v-i are the peculiar velocity 
potentials of CDM and DE, respectively, with the same order of 8\ and S 2 ; 8p^(t, x) and 
<5p^ ff (t,x) are the perturbations of x) and p^(t, x), respectively, 

^m( r > x ) = -o-[prf6(7-)<5 2 (r,x) - p mft (r)5i(r,x)], (61) 
£p^f (r,x) = <5p d (r,x) + o-[p d6 (r)5 2 (T,x) - p m6 (r)5i(r, x)], (62) 

where 5p d (r, x) = p rf (r, x) - p^r). 
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In the conformal Newtonian gauge, the first-order perturbed Einstein equations give 
us 1331 



?>Utf + k 2 ^ + m 2 (p = -4nGa 2 (5 lPmb + 5 2 p db ), (63) 
fcV + k 2 U<P = 47tGa 2 [( Pmb + p^)e x + ( Pdb + p£)0 2 ] (64) 

V>" + H(2i// + 0') + (2— - n 2 )<P + ^-(tf> - <f>) = 47rGa 2 (5p d (65) 

a 3 

V> - = 0. (66) 

Hereafter, primes denote the derivatives with respective to the conformal time r, TL' = 
a' /a, 0\ = —k 2 v\ and 9 2 = —k 2 v 2 . The first-order equation of the conservation law of 
CDM T^ e . b = (in Fourier space) tells us [23] 



s[ - -m^h + m— - 3(1 + ^)^' + (i + — )0i = o, (67) 

Pmb Pm?> Pmft Pm?> 

g/ + P'mb + Kg gi + _ k 2 (j) _ k 2 5Pm ^ = q (6g) 

Pmb + P m ft Pmb + P m fe 



The first-order equation of the conservation law of DE T%f!. b = tells us that [33] 

5' 2 - m^5 2 + m-^- - 3(1 + + (1 + — )0 2 = 0, (69) 

Pdb Pbb Pdb Pdb 

Pdb + p e d t pdb + vf b 



By differentiating Eq. flBTjl with respect to r and using Eqs fl65|) . (1661) . fl68l) and fl63|) . 
taking the geometric optic limit, finally we can get 

8'{ = -k 2 -^. (71) 

Pmb 

Similarly, from Eq.f l6"9~j) . and using Eqs ([6"5"j) . f l6"6"j) . ( TT0|) and f[6"3"j) . taking the geometric 
optic limit, we can also get 

S'» = -k 2 ^. (72) 
Pdb 

By using Eqs.f lSTj) and fl62|) . we can rewrite the above two equations as 

d 2 ,2w 



Here the matrix .M is 

M = ( _, , 

-or C^ + (T 



^ = ( l-i .2 ) > ( 74 ) 
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where 

_ Pdb 2 _ dPd 
r — > c sd — x ■ 

Pmb OPd 

This matrix possesses two eigenvalues A + and A_, 

A± = gL±j^5±jg. (75) 

Clearly, A± are crucial for determining the stability of the interacting model. In fact, X± 
are the physical sound speeds in the interacting model. If A + > and A_ > 0, the model 
is stable. Otherwise, the model is instable. Here we assume a > and it can be easily 
checked that c 2 sd > indicates A± > and c 2 sd < indicates A_ < 0. So we know c 2 sd is 
crucial for determining the stability of the interacting model, and negative c 2 sd indicates 
the existence of instabilities. 
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